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Abstract 

We compute the non-trivial infrared c^-fixed point by means of 
an interpolation expansion in fixed dimension. The expansion is for- 
mulated for an infinitesimal momentum space renormalization group. 
We choose a coordinate representation for the fixed point interaction 
in derivative expansion, and compute its coordinates to high orders 
by means of computer algebra. We compute the series for the criti- 
cal exponent v up to order twenty five of interpolation expansion in 
this representation, and evaluate it using Pade , Borel-Pade , Borel- 
conformal-Pade , and Dlog-Pade resummation. The resummation 
returns 0.6262(13) as the value of v. 



1 Introduction 



Non-trivial fixed points are a highly challenging aspect of renormalization 
theory. Much of what is known about non-trivial fixed point is due to the 
e-expansion of Wilson and Fisher ||WF72| , [WK74|| , which is an interpolation 
from a critical dimension to the one of interest. In this paper we present 
another interpolation scheme where the dimension of the underlying (Eu- 
clidean) space-time is kept fixed. 

A prototype of a non-trivial fixed point is the infrared fixed point of mass- 
less 4 -theory in three dimensions, also called Wilson fixed point |[WK74 . 



We choose it as an example for our method. Although it has been inves- 
tigated by various other means, for instance by Monte Carlo simulations 
of the three dimensional Ising model close to criticality, hopping parameter 
expansion, field theoretic perturbation theory for its scaling limit, and nu- 
merical integration of renormalization group flows in a number of setups, our 
knowledge of it is far from satisfactory. Accurate data for its spectrum of 
anomalous dimensions is lacking, its functional form is largely unknown, in 
particular its locality properties, and its mathematical construction remains 
an outstanding difficult problem. We mention [ [ZJ89|| and references therein 
as a guide to the extensive literature. We mention further that all this has 
been accomplished to a very satisfactory status in the hierarchical approx- 
imation by Koch and Wittwer ||KW91|1 . Our interpolation is a brick in the 
analysis of the full model. 

As starting point we choose a functional differential equation from the in- 
finitesimal renormalization group of Wilson ||WK74|| . Specifically we choose 
a normal ordered and rescaled representation for the fixed point interaction, 
expressed in terms of a scalar field with non-anomalous scaling dimension. It 
contains a bilinear renormalization form. This bilinear form is continuously 
turned on with an auxiliary parameter such that zero gives a linear theory 
and one restores the full equation. The linear theory is arranged such that 
the 4 -interaction acquires the scaling dimension zero in it. We then expand 
the fixed point interaction into a power series in the interpolation parameter. 
This part is similar to the e-expansion. In order to perform the expansion to 
high orders on the computer, the interaction is written in a basis of interac- 
tions which includes a general two point interaction in derivative expansion 
together with local higher interactions. We compute both the fixed point 
interaction and the eigenvalue associated with a massive perturbation. The 
resulting power series are evaluated by means of Pade , Dlog, Borel-Pade , 
and Borel-conformal-Pade resummation. 

The interpolation idea applies also to other renormalization schemes. The 
infinitesimal renormalization group is a particularly convenient one because 
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it involves a minimal set of Feynman integrals. Interpolations in fixed di- 
mensions can also be formulated for discrete renormalization group transfor- 
mations both in continuum regularization and on the lattice, at the expense 
of dealing with general non-linear rather than quadratic equations. It is con- 
ceivable that our interpolation can be given a meaning beyond perturbation 
theory. 

The paper is organized as follows. In section two we explain the structure 
of our particular functional differential equation. It is taken from ||W96|| and is 
the Wilson equation ||WK74j| in a kind of interaction picture. In section three 
we present our interpolation scheme. It is compared to a naive interpolation 
which has only a trivial solution. We solve our equations to lowest order 
to explain their recursive treatment. In section four we discuss their form 
in a coordinate representation. The result is a set of algebraic recursion 
relations for the fixed point interaction. They involve a set of structure 
constants whose computation again involves certain Feynman integrals and 
multiplicities. We devote section five to this issue. In section six and seven 
the eigenvalue problem for the scaling fields of the non-trivial fixed point 
and their anomalous dimensions is treated along the same lines. We restrict 
our attention to a mass perturbation of the fixed point and the associated 
critical index v. In chapter eight and nine the resulting recursions are studied 
by means of computer algebra. We conclude with a brief discussion of our 
results on the value v. 



2 Renormalization group fixed point 



We consider a real scalar field on three dimensional Euclidean space. We 
use a momentum space renormalization group built from the decomposi- 
tion of a massless propagator v with exponential ultraviolet regulator. The 
renormalization group will be formulated in terms of an interaction V(4>) . 
Concerning the general background, we refer to the work of Wilson ||W71|| , 
Wilson and Kogut [ WK74 |, and also to Gallavotti [G85| . Our setup will be 
identical with that in ||W96 |. 

We study the non-trivial infrared fixed point in three dimensions as so- 
lution to the functional differential equation 



v <t>>TI )V ((/>) = (V((f>)\V ((/>)), 



which was derived in [|W9£ |. Eq. ([!]) is a normal ordered and rescaled 



variant of the infinitesimal renormalization group due to Wilson | WK74 



Its origin is a flow equation governing the behaviour of interactions upon 
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the infinitesimal change of a floating cutoff. Eq. ([!]) gives stationary flows 
modulo the rescaling of units. The left hand side of (|I|) is a generator of 
dilatations 




L=l 



L-*t f I ) , (2) 



acting on the interaction. The field <fi is here rescaled non-anomalously with 
its canonical dimension at the trivial fixed point. It should be distinguished 
from the scaling fields of the infrared fixed point which have non-zero anoma- 
lous dimensions. The right hand side of ([!]) is a bilinear renormalization group 
form 




(V(<t>) \V(</>)) = — , x— exp — , v—\) VifaWfc) 




• (3) 

■4>2=4> 

It can be visualized as a sum of contractions between two copies of the 
interaction. Each contraction is made of one hard propagator \ and any 
number of soft propagators v. The propagators are here given by 

X{p) = e p , v(p) = — (4) 



as m 



W96|| . Eq. ([!]) is the main dynamical equation in this investigation. 
Being a differential equation, it has to be supplied with further data to select 
a particular solution. In a rigorous theory in the sense of Glimm and Jaffe 
[|GJ87|1 , the infrared fixed point should come as a global Z 2 -symmetric solu- 
tion, where global refers to some criterion of finiteness. Our point of view 
in this approach will be more modest. An interaction V(<j)) will stand for a 
power series 

oo „ 

V ( ( t ) ) = J2 d3x l ' ' ' 4>(Xl) ■ ■ ■ 4>(x2n) V 2n (Xl, . . . , X 2n ) (5) 

n=l J 

in the field, with symmetric Euclidean invariant distributional kernels given 
by Fourier integrals 

V 2n ( Xl ,...,x 2n ) = / . . . e ^ 1+ - +M n) 



(2tt) 3 (2tt) 3 
(27rf5 {3) ( Pl + ■ ■ ■ + P 2n)V 2n ( Pl , . . . ,p 2n ) (6) 

of smooth momentum space kernels. I.e., we identify an interaction with its 
collection of momentum space kernels. The question of convergence of the 
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expansion (|5|) in powers of fields will not be addressed. It is conceivable that 
it could be tackled with a suitable norm on the collection of momentum space 
kernels as a whole. 

In the iterative approach to be defined below we will meet at finite order 
no more than polynomial expressions in the field. We will understand ([!]) 
as a system of differential equations for the momentum space kernels. Its 
explicit form can be looked up in |[W96|| . Boundary data is substituted for by 
the condition of regularity. Homogeneous functions give particular kernels, 
which correspond to scaling fields of the trivial fixed point. Expanding a 
kernel in powers of momentum derivatives, we can always express it in terms 
of such scaling fields. To distinguish them from the scaling fields of the non- 
trivial fixed point and also because we will use perturbation theory, we will 
speak of them as vertices. 



3 Interpolation parameter and expansion 

Our strategy to solve (P is to interpolate to a solvable situation. If the 
interpolation is smooth it can be performed by means of perturbation theory. 
A natural candidate is 

(vc/>, Aj V(</>, z) = z (V(<f>, z)\V{<f>, z)) (7) 

with an interpolation parameter z — . . . 1. It can be thought to turn on 
continuously the bilinear form, which is identified as the source of troubles. 
The interpolation (0) is inappropriate for the following reason, when the 
dimension parameter is fixed to three. Expand the interpolated interaction 
as a function of the interpolation parameter in a power series 

oo 

v(M = £* p v< p >fo). (8) 

r=0 

Unfortunately there is little hope that (H) has a finite radius of convergence 
both in the case of (|?|) and the interpolation fll~5|) considered below. To 
be cautious we will therefore view (§) as a formal power series and inter- 
pret all equations below in this sense. It will however be argued that non- 
perturbative information can be extracted by Borel resummation. In order 
to solve (0), the expansion (||) has to satisfy 

[t>4>, ^) ^ (r) (0) = £ {v {s \<PW {r - l - s \<P)) (9) 
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to every order r G N, with the understanding ^(0) = 0. In particular it 
requires the interaction to satisfy 



V^) ^ (O) (0) = O 



(10) 



to zeroth order. In other words, the zeroth order has to be a marginal scaling 
field of the trivial fixed point. In three dimensions we have two marginal 
scaling fields, a wave function term and a 6 -vertex, to be abbreviated as 

0^(0) = J d 3 x 0(x)(-A)0(x), 3 ,o(0) = / d 3 x 0(x) 6 . (11) 

We emphasize that vertices should be understood as momentum space ker- 
nels at zero momentum and their Taylor expansions. Each of them comes 
with a formal orthogonal projector V\ } \ and V^o, selecting the corresponding 
vertex from a general interaction (|5]). The zeroth order has to be a linear 
combination 

= vf§ o.M) + v$$ e>3,o(0). (12) 



The coupling constants are not determined by the zeroth order equation (|10|) . 
To first order, (J^j reads 

(v^ j^j vv>(4>) = (y (o) (0)i^ {o) (0)> • (is) 

But eqs. (0) and (|13"D together have only a trivial solution. The left hand 
side of ( O ) cannot contain the vertices flTT|) because the dilatation generator 
(X>0, jhas no marginal image. Therefore, (|i~3"D requires that 

V hl (V<W(0|V< o) (0) = P 3 ,o (V (o) (0)|^ (o) (0)) = 0. (14) 

Computing the bilinear form with two copies of ([12]) inevitably gives V} ^ = 
= 0. Eq. ([7D is thus an inappropriate interpolation and has to be given 

up. 

A way around the obstacle is to interpolate simultaneously the dimen- 
sionality of the theory. This is the strategy of the e-expansion of Wilson 
and Fisher |[WF72|| in a field theoretic setup. Another way is to interpolate 



the scaling dimension, remaining firmly in three dimensions. We choose this 
second route and replace (|7|) by 



V(<j>,z) = z {V(</>,z)\V(<i>,z))-zV(<i>,z). (15) 
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The power series expansion ([8]) solves ( PUD if the coefficients satisfy the system 
of differential equations 



r-1 



= 2 (y w (0)|y (r_1 - ) (0)) - ^ (r_1) (0) (is) 



s=0 



to all orders r £ N. To order zero, (|f6D requires now 



0°)(0) = i 



(17) 



in contrast to (|H]). The zeroth order interaction is now a scaling field with 
unit scaling dimension. In three dimensions we have only one candidate, the 



-vertex 



O2M = / d 3 x 0(x) 4 . 



The zeroth order interaction thus has to be proportional to fll8"|) . The pro- 
portionality factor is the 4 -coupling. It is not determined by the zeroth 
order equation (P7|). We conclude that 



^ (O) (0) = Wo ) 02,o(0) 



(19) 



This expansion proves to have indeed a non-trivial solution. To see this, 
consider the first order equation in (TTH) . It reads 



V** 



v^\4>) = (v^(4>)\v {0) W)-v^\ ( p) 



(20) 



Eq. ( |20| ) cannot have a 4 -vertex on its left hand side. Therefore it is 
required that 

V 2 J (V {0) mv {0) (<!>)) - ^ (O) (0)1 = 0. (21) 



Computing the bilinear form, this condition reads explicitely 

y 2 (0 ) {l44x^^(0)^° ) -l} = 0, 



(22) 



where * means convolution times (2ty) 3 . Besides the trivial solution V 2 \°q = 
it has a non-trivial solution 



(0) 



1 



(2tt 



,3/2 



2.0 



U4x*v(0) 72 



0.21874445 



(23) 
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The value of the 4 -coupling at any given order will in fact be determined 
by the equations at the next order, a feature of this particular interpolation 
expansion. To first order, the interaction can be split into 



vW(0) = W 1 o ) O 2 ,o(0)+P 2 J - o r (1) (0) 



(24) 



with = 1 — 7\o the projector on the formal orthogonal complement. Eq. 



( |20D defines a system of first order differential equations for the momentum 
space kernels therein. They have a unique integral in the space of smooth 
functions of momenta, see ||W96 1 . We denote this integral by 



P 2 >«(0) 



v<f>, 



1 



vUv^\<t>)\v {0 \<i>) 



(25) 



This iterative scheme carries on to every order of interpolation expansion. 
Consider eq. (0) at order r > 2. The first step of the iteration is to 
compute V^q -1 '' at order r — 1. Making use of (|23|), its value follows from 



V, 



(r-l) 



2.0 



O2M 



2P 2 ,o(V {O) (0)l^2 J :o^ ( ^ 1) (0))- 

r ^v 2 ,Jv^\<t>)\v {r - l - s) {<t>))- 



s=l 



Again one splits the order r interaction into 

yM(0) = ^; o ) O 2jO (0)+P 2 >M(0) 



(26) 



(27) 



and computes the formal orthogonal complement to the 4 -vertex by inte- 
grating the first order differential equations (|16|). The result can be written 

as 



V<f>,T7 



' \xn fl (v^mv {r - i ~ s \<p) 

U=o 



:28) 



Thereafter it is time to proceed to the equations at order r + 1. For this 
scheme to work as above it is important that the kernel of (t>4>, j^) — 1 be 
one dimensional. Otherwise we would have to compute further order r — 1 
data from the equations to order r. An example where this happens is the 
4 -trajectory in four dimensions ||W96| . 



Although being in principle doable, the computation of this scheme to 
very high orders of interpolation expansion is a tedious enterprise. The main 
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work is the computation of a wealth of Feynman kernels generated in the 
course of iteration. A low order analysis of this program will be presented 
elsewhere. In this paper we choose to evaluate the expansion to high orders 
for a sub-class of contributions in the iteration. For this purpose we refor- 
mulate the fixed point equation (^) into an algebraic system of equations 
for a set of coupling constants. We find it interesting by its own. It also 
allows to perform the expansion on a computer. 




4 Coordinate representation 

We choose a system of vertices Oi(4>) labelled by elements i of an index set 
I. The vertices will be required to be Z2-symmetric, Euclidean invariant, 
and linearly independent. They will also be required to be regular in the 
sense that they are given by smooth momentum space kernels. We choose 
the system such that the dilatation generator acts linearly on it through a 
scaling dimension matrix 

CW) = i;<W)Ej. (29) 

We restrict our attention to systems with the property that the scaling di- 
mension matrix is diagonalizable. Non-diagonalizable matrices will not be 
considered here. In this case, we can arrange the system to consist of eigen- 
vectors 

(vfr-j^Oitt) = <?&(<!>). (30) 

In other words, we take Oi(4>) to be a scaling field of the trivial fixed point 
with scaling dimension <j{. Recall that such vertices are given by homoge- 
neous momentum space kernels. Eq. ( p9~|) says that the system closes under 
the action of an infinitesimal dilatation. We also require it to close under the 
action of the bilinear renormalization group form. For any two vertices Oi(<f>) 
and Oj(4>) the bilinear form (|3|) will be assumed to be a linear combination 

{o l moM)) = H°^) F t ] ( 31 ) 

feel 

with a set of structure constants F^. The scaling dimensions and the 
structure constants comprise all the information needed in the following 
about the system of vertices. We remark that the structure constants are 
well defined through 

F*O k (<l>)=V k (OiWOM)) (32) 
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even when the system does not close under the bilinear form. In this case 
(|3l|) holds only up to an error term. Below we will indeed work with an 
approximation of this kind and argue that the error term is small. 

We define a coordinate representation for the interpolated interaction in 
terms of a given system of vertices as 

v(M = £<W)^(4 (33) 



The idea is then to investigate the interpolation fll5|) for the infrared fixed 
point by means of the parameter dependent coordinates (0). Eq. (|T5|) 
becomes a system of algebraic equations 

(a k - 1) V k {z) =zJ2 F t,j V*(z) V(z) - z V\z) (34) 
ijei 

in the coordinate representation. The advantage of (|34"D as compared to (|15f) 



is that we are no longer dealing with differential equations for momentum 
space kernels. Their integration is hidden in the structure constants. If the 
interpolation is smooth, we can expand the coordinate functions into power 
series 

oo 

V k (z) = £ z r yk ' r - (35) 

r=0 

By standard arguments Q3"5D is expected to be singular but Borel summable. 
Our below evaluation of fljffi) supports this expectation. Eq. (^5]) yields a 
solution to fl3"4|) in the sense of a formal power series in z if the coefficients 
obey 

r-l 

(a k - 1) V k ' r = Y.Y. F ti yi ' S V 3 ''"- 1 - 3 - V^- 1 (36) 

s=0 i,j(zl 

holds for all couplings k G X to all orders r 6 N of interpolation expan- 
sion. We organize (^) into a recursion relation which can be solved on the 
computer. To zeroth order (^) simplifies to the linear equation 

K - 1) V k '° = 0. (37) 

We assume that our system of vertices contains only one element labelled by 
k = 2 = (2, 0) such that oi — 1. This element is of course the 4 -vertex 



(|18). All other elements are assumed to have scaling dimensions different 



from one. Then ([J7]) implies that 
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The value of V-'° is as above determined by ( |36| ) to order one, 

i?k - 1) V k > 1 = (Fi >2 _ V 2 -' - 5 1>k ) . (39) 



Evaluating (|39|) for k — 2 it follows immediately that we have 

V 2 -' = 4-, (40) 

besides the trivial solution V^-' = 0. Eq. (^) does not tell the value of V-' 1 . 
But for k G X \ {2} it gives 

,1 *2 2 (^-'°) 2 

CTjfc - 1 

Eq. ([40]) and (|41~|) are of course the coordinate versions of (|23|) and (p5|) . The 
strategy to any order r > 1 is again to first compute V-' r ~ l and thereafter 
^fe,r f Qr k £j\ {2}. The explicit formulas are 

r-2 

7 2,r-l = _ 2 ^ ^/2,0 yi.r-1 - £ £ V*'* 7 i ' r_1 -* (42) 

*GI\{2} ' s=li,jeT 

and 

1 f r 1 1 

= ) y~] F k - v t,s v j,r ^ 1 ^ s — V k,r ~ 1 > (43) 

in complete analogy to and (E3). Thus once we know the scaling di- 



mensions and the structure constants, the iteration proceeds by means of 
purely algebraic operations. We remark that the sums in (|42"D and ([43|) will 
be finite in the system of vertices considered below. The reason is that the 
outcome of the bilinear form of two monomials in the field is a polynomial 
in the field of finite order, and consists only of connected vertices. A very 
interesting question is whether it is possible to find finite systems of vertices 
that close under both ([29]) and ([31]). It is clear that this cannot be achieved 
in terms of polynomial vertices. Unfortunately no such system is known in 
three dimensions. 

5 Structure constants 

We consider the following system of vertices. First we include a full two point 
vertex in derivative expansion. A convenient notation for it is 

(9i, a (0)= /d 3 x0(x) (-A) Q 0(x), (44) 
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where a = 0,1, 2,.... Second we include local (2n)-point vertices with arbi- 
trary many external legs. They will be abbreviated as 



a,o(0) = / d 6 x (j>(xY n , (45) 



where n = 2,3,4,.... Notice that both (44) and (|45| ) meet the demands 
stated at the beginning of the previous section. More general interactions 
include also momentum dependent higher vertices (|45|) . They will not be 
considered here. Our index set is thus 

J = {1} x {a e N\a >0}U{nG N|n > 2} x {0} (46) 

and 2 = (2,0). The bilinear form does not close under this set of vertices. 
For instance two local vertices fl4"5| ) contract in general to a bilocal vertex. 
Thus if we perform an iteration (|42| ) and (|43"D with this system of vertices, 
we make a systematic error due to the truncation of the system. Our ansatz 
rests upon the assumption that non-local higher vertices are small compared 
to their local parts. 

The scaling dimensions of ( (44| ) and (|45|) come out as 



a 1>a = 2 - 2a, a nfl = 3 - n. (47) 

The structure constants for this set of vertices come out as follows. Two 
quadratic vertices always contract again to a quadratic vertex. The associ- 
ated structure constants are computed to 

^ ^ 1 \ a— 3 —~: 

F (iMi,<r) = 4 ( a _£_ 7 ), e ^+7' (48) 

where = 1 for a > b and zero else. A quadratic vertex and a higher vertex 
return upon pairing both a quadratic vertex and a higher vertex. First we 
have 

4 1 ;5,( 2 ,o)=24^i, Q (0)5 7 , . (49) 
The structure constant (^9|) involves the one loop integral 



KiJp) = / (^(P - ?) ( 5 °) 

at zero momentum. Recall that the propagators are given by (f|). The 
exponential regulator gives a convergent integral which is evaluated in fl5"8"p. 
Second we have a one loop contribution 

^fcjSo) = 4 ™( 2 ™ - !) ^ (°) ( 51 ) 
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as well as a zero loop contribution 

^ko) = 4m <W" ( 52 ) 

This last pairing also contributes to momentum dependent higher vertices 
which we neglect. Two higher vertices yield upon pairing both a quadratic 
vertex and higher vertices. One quadratic term is 

4^,0) = - l)(2n)! Sg_ 2 (0) (53) 

with the (2n — 2)-loop (the number of soft propagators v) integral 

K 2 „-2(p) =v*---*v*x{v) (54) 

expanded into 

oo 

K 2n . 2 (p) = Y,(p 2 TK^ 2 (0) (55) 
at zero momentum. A second quadratic term is 

<o),(n-i,o) = (n - 1)(2")« *2n- 2 (0) V (56) 

This second term is exactly local. The general higher vertex content of the 
pairing of two higher vertices is summarized in 

F (i,o) (2n)! (2m)! 

(n,0),(m,0) ( n + m _/_ l)!( m + Z- n )!(Z + 

n — m ! 

-^n+m-Z-l(0) O n + mj i+i O m +/ jn G; +rljm . (57) 



This last set of structure constants ( |57| ) alone defines a local approximation 
for the renormalization group fixed point. As mentioned above, the gen- 
eral outcome of the pairing of two higher vertices also contains momentum 
dependent vertices which are not encorporated in fl5"7|) . All other structure 
constants between vertices in X are zero. 
The one loop integral fl5"0|) is evaluated to 



^i, a (0) = -^2 1/2 - a r(a + l/2). (K 



The /-loop integral fl5"j|) is computed as a function of the external momentum 
squared to 

% (p) = (4tt) - 3 '/ 2 dai • • • da, A" 3/2 exp (^V ) (59) 
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with the abbreviations 

1+1 i 
A = E II a n, B=l[a m , (60) 

m=l n^m m=l 

where ai + \ = 1. Its momentum derivatives at zero can be reduced further to 
a one dimensional integral 



K/ (U) - (Stt)' a!I> + 3/2) Jo I J 



(61) 



This remaining integral can be done explicitely at least in the one-loop case. 
We evaluated it in the general case numerically to high accuracy (45 digits) 
on the computer. 



6 Eigenvalue problem for critical indices 

The fixed point equation ([[]) comes together with an eigenvalue problem 



V(f>, 



A 



W(<f>)=2 (V{cf))\W(<f>)) 



(62) 



defining scaling fields W(4>) and their anomalous dimensions A. We em- 
phasize that W(4>) is a composite field of (f>. The spectrum of anomalous 
dimensions is the object of principle interest associated with a fixed point. 
It directly determines the critical exponents, see Wilson and Kogut jWK.74 . 
The interpolation ( |i5| ) is accompanied by an interpolation of fl6"2|), given 

by 

* 1 • ' W&z) = 2z (V(<P,z)\W^,z)). (63) 



X(z) 



Eq. ( |6"3"D can be solved by means of perturbation theory. We expand not only 
the interaction (|8|), but also the scaling field and its anomalous dimension 
into power series in the interpolation parameter, 



(64) 



r=0 



r=0 



We interpret ( p4|) in the sense of a formal power series. It yields a solution 
to fl63|) if the coefficients satisfy the system of differential equations 



S 



r-l 



V(j>,j-\ W {r \ct)) -^A {s) I^ (r - s) (0) = 2^(y (s) (0)|^ (r_1 " s) (0))- (65) 



s=0 



s=0 
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This system can be integrated iteratively. To order zero, (|65| ) becomes the 
eigenvalue problem 




A (o) 



W {0 \4>) = 0. (66) 



The zeroth order W^((p) thus has to be a scaling field of the trivial fixed 
point, and A 1 - ** has to be its scaling dimension. With each scaling field of 
the trivial fixed point is therefore associated in perturbation theory a scaling 
field of the non-trivial fixed point. 

Let us consider for definiteness the perturbation associated with a mass 
term 

W^{<f>) = llO (0), lfi {4>) = I d 3 x 0(x) 2 . (67) 



Then the zeroth order eigenvalue is of course 

A (o) 

= trio = 2. As a per- 



turbation of the non-trivial fixed point, (|67|) turns out to be relevant. The 
associated non-trivial renormalized trajectory in the sense of ||W96|| describes 
the renormalization group flow of a non-trivial massive field theory. Associ- 
ated with it is the critical exponent 

v = i. (68) 

The mass perturbation flBTD is non-degenerate in the sense that the kernel 
of (v<j>, ±) - 2 is one dimensional. The formal orthogonal projector on this 
one dimensional kernel is V\$- Another non-degenerate perturbation is the 
scaling field associated with C 2iO (0). The ones associated with 0\^{§) and 
C^3,o(0) on the other hand form a degenerate marginal duplet. We will restrict 
our attention to the non-degenerate case for the sake of notational economy. 
The kernel of (T>4>, j^j — will thus be assumed to be one dimensional. 
The formal orthogonal projector on this rank one kernel will be denoted by 
V . To first order (|65[) becomes the differential equation 

W w ((f)) - \ w W i0 \<f)) = 2 (v^((j))\W^((j))) . (69) 




A (o) 



The first order correction to the eigenvalue follows from (|69|) by projection 
with V. We have that 

\V) W (o) (0) = -2V (V (o) (0)|^ (o) (0)) • (70) 

Spelled out explicitely for the mass perturbation (]67f) , eq. (|70|) says that 

\V = -4&x*mV$l = I £- (71) 
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It is amusing that this first order correction can be inferred without having 
to compute the convolution integral, because the convolution integral in ([711) 
is canceled exactly by the one in (E^). Next we impose the normalization 
condition that 

VW {l \<p) = Q. (72) 

This condition is appropriate in the non-degenerate case because ( |70D al- 
ready takes care of the P-information contained in fl6"9|). The orthogonal 
complement is then integrated as in the case of the interaction. The out- 
come is 



V L W {x \(j)) 




_ A (o) 



-i 



v ± (v (0 \4>)\w^(4>) 



(73) 



This scheme carries on immediately to every order of interpolation expansion. 
Projecting (^) to V, we first deduce that 



r-l 



A (r) W {0) (0) = 2 J2 V (V (s) (0) | W^- 1 -) (<p) 



s=0 



(74) 



This equation determines the order r eigenvalue in terms of lower order data. 
Generalizing (|72|), we impose the normalization condition 



VW ir \(f>) = 

for r > 1. Then to order r we are left with the computation of 

5 



(75) 



V<f>, 



r-l 



A (0) 



-1 (r-l 



Y / \ {s) v ± w (r - 1 \0) + 



.s=l 



2 E^ (^ (s) (0)l^ (r-1-s) (0))- 



s=0 



(76) 



This scheme iterates to every order of interpolation expansion. Recall that 
the inverse of the dilatation generator in (|76|) involves the integration of a 
first order partial differential equation. As in the case of the fixed point, the 
explicit computation of this program to very high orders requires considerable 
computational resources. In this paper we restrict our attention to a partial 
resummation by means of our coordinate representation. 



7 Eigenvalue problem in coordinates 

In this section we perform the interpolation expansion for the eigenvalue 
problem ( |5"3D in the coordinate representation. The coordinate representation 
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for the scaling fields reads 

w(M = £<W)^(*)- ( 77 ) 

In the coordinate representation, the eigenvalue problem ( |63"D becomes a set 
of algebraic equations 

K - A(z)) W k (z) = 2zJ2 F t k ;j V\z) W*(z). (78) 

ijex 

It can be solved recursively in an interpolation expansion 

oo 

w k (z) = J2* r wk,r - ( 79 ) 

r=0 

The power series fl35|), (|64|), and ( |79|) yield a solution to ( |7H|) provided that 
the coefficients satisfy 

r— 1 i — 1 

(cr fc _ v^ fc ' r - A (r) w K0 = w Kr ~ s + F ij yi,s W^" 1 ' 8 . 

s=l s=0i,jel 

(80) 

As in the case of the interaction, the system of equations ( |8"0"D can be orga- 
nized into a recursion relation. To order zero (^) reads 

(a k - A (0) ) W k '° = 0. (81) 

It tells us that we should select one of the k e I as zeroth order eigenvector. 
We choose 1 = (1, 0) for defmiteness. Then the zeroth order is 

W k '° = S Lk , A (0) = (7i = 2. (82) 

The only k with a k = 2 is k = 1. We will again restrict our attention to this 
non-degenerate case. The below recursion relation is valid for general non- 
degenerate perturbations, with minor notational changes. The first order 
equation in the system ( |8~0"D is given by 

(a k - A<°>) W k,x - \Ww k '° = 2 F^V^. (83) 

Therefrom it follows that the first order correction to the eigenvalue is in the 
coordinate representation 

A« = -2 F^V g '°. (84) 
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We remark that in the degenerate case, the degeneracy is typically lifted by 
the first order correction to the eigenvalue. The other coefficients to first 
order are 

W 1 ' 1 = (85) 



and, for k el \ {!}, 



rkA Z ^2,1 



This computation generalizes immediately to higher orders. The formula for 
the order r eigenvalue in terms of lower order data is 



\(r) = -2j2Yl F h yl ' S W^- 1 - 3 . (87) 

s=0 



The order r eigenvector is then given by 

W l ' r = 

for r > 1, together with 

'r-1 r-1 

y a (s) iy fc - r - s + 2 y 



H/fe ' r = ~ _A(0) 1 S aW W/fc ' r ~" + 2 EE^ ^' r - 1 - s [, (89) 



for k G X \ {1}. Eq. (|87D , (|88|) , and (|^) define a recursive perturbation 
expansion for the critical indices of the non-trivial fixed point. 



8 Computation of the recursions 



We computed the ^-expansion for the potential recursively by means of (|42|) 



and (pESp, and for the eigenvalue problem by means of ([37]) and (^) using 
computer algebra. We restricted our attention to the case of three dimen- 
sions. It turned out to be crucial to compute the structure coefficients to 
high accuracy. We calculated them to an accuracy of 45 digits with Maple 
V. The perturbation expansion was performed up to a maximal order of 25. 
The derivative expansion was performed up to = 20 orders of p 2 in 
the 2-point vertex. Table |l] shows the series for the 4 -coupling both in 
the ultra-local approximation = and for = 4 up to the order 
z 11 . The coefficients prove to increase in absolute value proportional to C n n\ 
with some constant C . Their signs alternate. From this behavior we conclude 
that the series does not converge but is Borel summable. A proof of local 
Borel summability will be presented elsewhere. The constant C is related to 
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n 




V ( } cv — 4 

^2,0 5 "mat ~~ ^ 





2.1874-10" 1 


2.1874-10" 1 


1 


4.5814-10- 1 


4.5814-10- 1 


2 


-8.7171-10- 1 


-8.6761-10- 1 


3 


4.6575-10° 


4.6815-10° 


4 


-4.0553-10 1 


-4.0546-10 1 


5 


4.2980-10 2 


4.2992-10 2 


6 


-5.2117-10 3 


-5.2130-10 3 


7 


7.0118-10 4 


7.0133-10 4 


8 


-1.0267-10 6 


-1.0269-10 6 


9 


1.6155-10 7 


1.6158-10 7 


10 


-2.7080-10 8 


-2.7084-10 8 


11 


4.8059-10 9 


4.8066-10 9 



Table 1: Examples for the behaviour of the expansion coefficients. 



an instanton singularity of the Borel transform on the negative real axis. It 
can be seen as an accumulation point of poles when the series is converted 
into various Pade approximants. The derivative expansion on the other hand 
proves to converge. This is illustrated in table || for two values of aw- We 
note in passing that the difference between aw = 5 and a,^ = 10 is small. 

The spectrum of the non-trivial fixed is computed along the strategy ex- 
plained in section ^ and [7]. It requires as an input the fixed point interaction 
in z-expansion. We evaluated it for all values of am^ inbetween zero and 
twenty. In the following we will concentrate on an estimate of the critical 
index v (|68f) by resummation of the series for all these twenty one approxi- 
mations. We computed the series by means of (|57p and ( |59"D to order twenty 
five of z-expansion. 

Table ^ shows as an example the series for the eigenvalue A in the ultra- 
local = and in the case of a^ = 4 up to the order twelve of 
perturbation theory. Again the series alternate, and the coefficients grow in 
absolute value as C n n\. The series are therefore not expected to converge. 
We remark that a proof thereof is however missing. The Borel transform 
of a series with this asymptotics has a finite radius of analyticity R armx - It 
is determined by an instanton singularity on the negative real axes of the 
complex Borel plane. This radius of analyticity is an interesting quantity. 



It can be investigated by a number of methods, see ||DI89|| and references 



therein. One of them is the Pade method. Recall that the Pade approximant 
of order (/, m) for a function / is a rational function fi jm (z) = ■ Here Pi 
and Q m are polynomials of degree / and m respectively, determined such that 
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a 


T/ ( 0) rv - ^ 
v 2,a ) "mar — <J 




= 10 





-1.08730107-10 5 


-1.08730280 


10 5 


1 


-7.30254527-10 4 


-7.30254464 


10 4 


2 


9.42673139-10 3 


9.42673053 


10 3 


3 


-1.47739400-10 3 


-1.47739385 


10 3 


4 


2.14057814-10 2 


2.14057791 


10 2 


5 


-2.79810089-10 1 


-2.79810056 


10 1 


6 


0.00000000-10° 


3.33201984 


10° 


7 


0.00000000-10° 


-3.70314164 


io- 1 


8 


0.00000000-10° 


3.99772841 


io~ 2 


9 


0.00000000-10° 


-4.43094961 


10~ 3 


10 


0.00000000-10° 


5.30798357 


IO" 4 



Table 2: Examples for the behaviour of the mass coefficients at order z 
for dmo,. = 5 and a m „ = 10. 



n 


\ n rv — n 


A -i Ctmox 


= 4 





2.000000-10° 


2.000000 


10° 


1 


-3.333333-10- 1 


-3.333333 


io- 1 


2 


-3.490659-10- 1 


-3.490659 


IO" 1 


3 


1.148993-10° 


1.159189 


10° 


4 


-7.414413-10° 


-7.369227 


10° 


5 


6.358855-10 1 


6.358630 


10 1 


6 


-6.649232-10 2 


-6.646081 


10 2 


7 


7.999490-10 3 


7.996744 


10 3 


8 


-1.070838-10 5 


-1.070532 


10 5 


9 


1.562548-10 6 


1.562166 


10 6 


10 


-2.452524-10 7 


-2.452002 


10 7 


11 


4.103373-10 8 


4.102601 


10 s 


12 


-7.271917-10 9 


-7.270698 


10 9 



Table 3: Series coefficients for A up to order 12 for = and aw = 4. 
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the taylor expansions of / and fi iTn agree up to order z +m . One then observes 
that the poles of the various possible Pade approximants accumulate around 
a cut or a singularity of /. With the Pade method we found 

R a _= 0.88 ±0.02, (90) 

with no significant dependence of aw. Figure [I] shows a plot of all poles of 
all Pade approximants (BX)i }jn with l+m = 25 in the complex Borel plane for 
the two cases aw = and a^ = 4 respectively. Here B\ denotes the Borel 
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Figure 1: Radius of convergence of BX by the Pade method for a^ = and 



transform of A. As expected, the poles accumulate on the negative real axes. 
Notice however that there are many spurious singularities on and nearby the 
positive real axes. These spurious poles endanger the inverse Borel transform 
as a contour integral along the positive real axis. The pictures for a^, = 
and oimuc = 4 show tiny differences. For instance, the poles on the positive 
real axis are not on fixed locations and can therefore be regarded as spurious. 



9 Determination of v 



To compute the value of the critical index u, we have to evaluate the z- 
expansion at z = 1. Naive evaluation does not give a meaningful answer since 
the expansion does not converge. Therefore we had to rely on resummation 
technology. A review of series resummation and references to the original 
literature is given in ||ZJ89|| and ||DI89|| . We tried four standard methods and 
compared the results. 
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First we computed (for all values of a niax between zero and twenty) all 
Pade approximants (X)i, m with / + m < 25, and evaluated them at z = 1. 
These values are conveniently displayed in a Pade table (l-m grid). To get 
an idea for the value of A and an estimate for the error we computed the 
mean value and deviation for the lines of fixed order in z (I + m = const) 
in these diagrams after having discarded all values below a lower value \ m in 
and above an upper value \ max - The idea thereof is that large deviations 
come from spurious singularities. We were careful not to choose the window 
too narrow. Our error estimate should be regarded as rather pessimistic. If 
these mean values converge with increasing order in z we use them and an 
inspection of the whole table to find an estimate for the value of v. 

In the second method (Dlog) one computes the Pade approximants for 
the logarithmic derivative -jhj- A is then reconstructed as the exponential 
of an integral 

A^ = A(0)e^ 2 ^ logA( <™. (91) 

The integration can be performed numerically to high accuracy. The Dlog 
method is particularly efficient when the singularity is of the type X(z) = 
^ x _^ x ^ with a nonintegral exponent 7. 

The third proposal is to use a Pade approximants for the Borel transform 
of the series. The Borel transform of a power series f(x) = Y, n >o fn% n is 
defined by 

(Bf)(*) = T,&- ( 92 ) 

The Borel transform of power series with finite radius of convergence defines 
an analytic continuation of the function to a maximal simplex through the 
integral 

roo 

f(x) = / dte-\Bf)(xt). (93) 







Again we get Pade tables of approximants for A by numerically integrating 
this back transformation for various Pade approximants of BX. 

The off diagonal estimates in these tables can be improved by using in- 
formation on the analyticity properties of the Borel transform. Bf could for 
instance have a cut along (—00, —R] on the negative real axes. Let us assume 
that this is indeed the case (with R = 0.88 ± 0.02). Then the cut plane can 
be mapped conformally via 

Jz/R+1 - 1 ARu 
Jz/R +1 + 1 (1-u) 2 
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onto the unit circle. Under this mapping (Bf)(z) transforms to (Bf)(u). 
We then use Pade approximants for the mapped series. A Pade table for A 
is obtained via the inverse transformation 



i 



l + u 




AR 



du 



(l-«) 



(Bf)l,m(u) . 



(95) 



The outcome of this method relies on a careful estimate of the radius of con- 
vergence of the Borel transform. For each am^, we calculated three estimates 
for A, one for our estimated value of R and one for R + AR and R — AR 
respectively, where AR means the error in our estimate for the error of the 
radius. The inspection of all three Pade tables yields A and an error estimate. 

We also tried out inhomogenous differential approximants, but we could 
see no improvement as compared with Pade or Dlog Pade approximants. The 
integration of the differential equations in this method turned out to be both 
time consuming and fragile due to the poles close to the origin. 

In table § we summarize our results for v for the different values of a^j. 
The errors refer as usually to the last digit. We come to the following con- 
clusions. The Pade method is the least precise one with an error of about 
0.01. From it we can get an idea about the value of u, but no accurate 
estimate. The errors of the Dlog Pade method (Dlog) and the Borel Pade 
method (BP) are of comparable size. The Borel Pade method with conformal 
mapping (BPconf) has the least errorbars. At higher orders of the derivative 
expansion of the 2-point vertex the errors increase significantly. To display 
this effect, we have plotted the data of table f| in figure ||. One can see that 
the values for v oscillate around a mean value. Up to a certain order this 
sequence seems to converge. Thereafter, the difference between u(a rmc ) and 
v{a-mac + 1) and the error grows. 

We believe this effect to be the consequence of a numerical instability. In 
high orders of derivative expansion and high orders of perturbation theory 
one is dealing with numbers of enormously varying magnitudes (in our case 
one hundred orders). In practice we computed our series to an accuracy 
of 45 digits, and a problem arises in the cancellation of large numbers in 
the course of the recursion. A more destructive explanation would be that 
the resummation fails to produce a convergent derivative expansion, or even 
more desastrous that the non-perturbative kernels are not analytic functions 
of the momenta. The final answer to this question can only be given on the 
basis of a non-perturbative construction of the fixed point and is outside 
the scope of this paper. Our insight comes from the evaluation of various 
approximants to different orders of accuracy. 

We confine our further discussion to those values of which lie before 
the onset of instability. The Dlog method and the BPconf method both yield 
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v Pade 


7 / 111 t~\ 

V UlOg 


V Dr 


v BPconf 


n 
U 


U.oooU^zU J 


U.004U^1U ) 


u.ooyy^ou ) 


U.OOoU^oU J 


1 


U.OzUU^lOU J 


U.oloU^OU ) 


U.OIOU^ 1 v) 


U.01UU(1U ) 


o 

z 


n fiQ/i rif'i Qd^ 


U.ozyU(4U J 


U.ozo4^1o J 


U.ooUU^lU ) 


Q 

o 


U.ooUU^llU J 


U.ozzU(4U J 


n coon/' 7rA 


U.ozUU(4U J 


/I 
4 




U.ozoU^oU J 


U.ozoo^Oo J 


U.ozoU^lU ) 





U.OOUU^lUU J 


U.OZOU^4U ) 


U.OZ4U^OU ) 


U.OZOO^lU ) 


6 


0.6330(100) 


0.6260(40) 


0.6290(40) 


0.6270(10) 


7 


0.6310(110) 


0.6280(70) 


0.6220(60) 


0.6260(20) 


8 


0.6330(90) 


0.6230(70) 


0.6310(60) 


0.6266(4) 


9 


0.6290(90) 


0.6260(80) 


0.6230(60) 


0.6220(40) 


10 


0.6330(90) 


0.6260(40) 


0.6320(70) 


0.6286(3) 


11 


0.6300(110) 


0.6260(50) 


0.6200(150) 


0.6310(50) 


12 


0.6310(150) 


0.6260(60) 


0.6340(60) 


0.6305(25) 


13 


0.6280(120) 


0.6260(80) 


0.6200(200) 


0.6350(60) 


14 


0.6360(110) 


0.6270(50) 


0.6410(80) 


0.6440(30) 


15 


0.6330(100) 


0.6250(60) 


0.6200(180) 


0.6440(30) 


16 


0.6340(200) 


0.6310(40) 


0.6520(50) 


0.6300(50) 


17 


0.6250(150) 


0.6270(200) 


0.6320(200) 


0.6259(40) 


18 


0.6380(60) 


0.6380(110) 


0.6549(58) 


0.6590(160) 


19 


0.6420(150) 


0.6000(400) 


0.6060(130) 


0.6240(70) 


20 


0.6420(80) 


0.6420(200) 


0.6550(40) 


0.6430(160) 



Table 4: Results for the critical exponent v with the Pade method, the 
Dlog Pade method (Dlog), the Borel Pade method (BP) and the Borel Pade 
method with conformal mapping (BPconf) for various orders of the derivative 
expansion of the 2-point vertex. 
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Figure 2: Tie critical index v as a function of the order of the derivative 
expansion of the 2-point function for the four extrapolation methods. 
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nearly constant values for v at orders between a r 



4 and a r 



12 and 



between = 4 and aw = 8 respectively. We propose this value to be the 
limit of v at arbitrary order of the derivative expansion. 

Consider the data for the ultra-local case aw = and to first order 
a max — 1 of derivative expansion. For the ultra-local case, which can be 
compared with the hierarchical model [y = 0.6501625, | KW88|| ) we find 
v = 0.6625(33) which is bigger than the full critical index. Disregarding the 
pure Pade estimate, we get for a m = 1 the result v = 0.6144(62). This 
value is considerably lower than the value at a,^ = and even lower than 
the full critical index. In view of the tiny differences between the fixed point 
coefficients at a„ = and a m ^ 0, we find this surprising. Compare for 
example the coefficients in table III. 

At higher orders of derivative expansion, the values for v oscillate and 
converge to a mean value. The limit value has been determined as the mean 
values of v over the nearly constant plateaus. As best estimate for the BP- 
conf method we get v = 0.6262(13). The Dlog method yields v = 0.6259(57). 
These results should be compared with the critical index v of the three di- 
mensional Ising model in the literature. In table [5] we list a few results for v. 
A comprehensive article on this issue is ||BLH95|| . It also contains an overview 



V 


Method 


Literature 


0.6300(15) 

0.6298(7) 

0.630 


three dimensional renormalization group 


[ 


i ;z.n 

BB85 

|N91| 




0.6305(25) 


renormalization group, e-expansion 


1PZJ85| 


0.6301 
0.6300(15) 


high temperature series 

high temperature series for bcc-grid 


|K95| 
[NR90 


1 


0.6289(8) 
0.6301(8) 


Monte-Carlo methods 


[ 


FL91 
BLH9 


1 


0.625(1) 


Monte-Carlo renormalization group 


IMT9SII 


0.626(9) 


Scaling-field method 


|NR84j 



Table 5: Results for the critical exponent v of the full model. 



of experimental data. With series expansion and Monte-Carlo methods one 
gets v = 0.630. On the other hand the Monte-Carlo renormalization group 
suggests v = 0.625. This gap is object of current discussions. Our value is 
closest to the value of ||NR84|| and ||CT96|| . 
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10 Summary and discussion 



In this article we investigated a form of Wilsons infinitesimal renormalization 
group. The starting point was equation ([!]). We found a practical way to solve 
the equation in a systematic manner. The central idea was to introduce an 
interpolating parameter z, which continuously turns on the non-linear term 
in (|I|). Everything was expanded in this parameter. The interpolation was 
arranged such that the zeroth order is a 4 -vertex. The expansion was pre- 
sented both in a coordinate free representation and in coordinate form, where 
the interaction is expanded in a basis of vertices. As a basis we advocated 
the use of a full two point interaction in derivative expansion together with 
local vertices of any power of fields. Derivative interactions of higher powers 
were neglected. The basis of interactions came encoded in a system of scaling 
dimensions and structure constants. Their evaluation was reduced to a one 
dimensional Feynman integral which we evaluated numerically. We reformu- 
lated our expansion into recursive equations for the fixed point interaction, 
its scaling fields, and their anomalous dimensions. We performed a detailed 
analysis of the series for the critical exponent associated with a massive per- 
turbation of the fixed point. The result is a new and independent calculation 
of the critical index v of the the three dimensional Ising model. We solved 
the recursion relations for the eigenvalue problem up to high orders and an- 
alyzed the resulting series by means of four different extrapolation methods. 
Our best estimator for the critical index v is v = 0.6262(13). We compared 
our results with values for the critical exponent v known in the literature. 

The results encourage us to further investigations. On the menu of open 
problems we have the inclusion of momentum dependent higher vertices for 
the scalar model, theoretical estimates on the z-expansion, and the general- 
ization to vector and matrix models. We hope to return with accurate data 
on their critical properties by means of z-expansion in the near future. 
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